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Abstract 

We compute the populations of isotopes of Boron, Carbon and Nitrogen 

measured experimentally in intermediate energy heavy ion collisions. A two 

component soluble statistical model is used to find the initial populations of 

different nuclei at a finite temperature. These initial populations are both in 

particle stable and particle unstable states. The particle unstable states then 

decay. The final populations after these decays are computed and compared 

with experimental data. 
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I. INTRODUCTION 



In this work we attempt to calculate the populations of various isotopes of Boron, Carbon 
and Nitrogen that were measured in a number of experiments at the NSCL-MSU facility [1- 
2]. The calculation proceeds in two stages. In the first, primary populations are calculated 
in a two-component statistical model. The calculations in the first part are exact although 
numerical. These populations are both in particle stable and unstable states. In the second 
stage the particle unstable states are allowed to decay. This is done in a Weisskopf formalism. 
Exact calculations are very long and some approximations had to be introduced. These 
approximations will be discussed. After the decays, the populations are compared with 
experiments. 

One motivation for this calculation was that it serves as an application of the two- 
component statistical model where an exact calculation can be done. This, therefore, could 
serve as a benchmark of how far one can trust the predictions of the model. Unfortunately, 
the particular predictions we are looking for are also affected by the subsequent decays. This 
effect is not small. Hence, the predictions are the result of the combination of two models 
which had to be applied in tandem before experimental data could be compared. A recent 
application of the two-component model was the computation of the caloric curve || in 
nuclei. 

The sections are organised as follows. Section II gives a brief description of the two- 
component statistical model. After presenting, in section III, in words and simple formulae, 
the overview of the secondary decay calculation, we present in section IV the formalism that 
we use to model secondary decay. In section V we present some calculational details, section 
VI presents the results of the calculation. Summary and conclusion are presented in section 
VII. A short appendix of the more complicated formulae are presented in section VIII. 
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II. THE TWO COMPONENT SOLUBLE STATISTICAL MODEL. 



For completeness, we present here the essential details of the two-component statistical 
model. The one component model was described elsewhere The formalism of the 

two-component model can also be found in H 

Assume that the system which breaks up after two heavy ions hit each other can be 
desribed as a hot, equilibrated nuclear system characterised by a temperature T and a freeze- 
out volume V within which there are A nucleons (A = Z + N). The partition function of 
the system is given by 

Qz,n Xl"-,^T (2- 1 ) 

i<>i,j- 

Here n it j is the number of composites with proton number i and neutron number j, 
and u>ij is the partition function of a single composite with proton, neutron numbers i,j 
respectively. There are two constraints: i n i,j — Z an d J n i,j — N. These constraints 
would appear to make the computation of Qz,n prohibitively difficult, but a recursion rela- 
tion exists |3|,[J which allows numerical computation of Qz,n quite easy even for large Z and 
N. Three equivalent recursion relations exist, any one of which could be used. For example, 
one such relation is 

Qz, n = - i^i,jQz-i,n-j (2-2) 
Z id 

All nuclear properties are contained in cjij. It is given by 

^• = ^(^) 3/2 (^+J) 3/2 x^ (2.3) 

Here Vf is the free volume within which the particles move; Vf is related to V through 
Vf = V — V ex where V ex is the excluded volume due to finite sizes of composites. We take 
Vf to be a variable of the calculation, it is set to be equal to fVo where Vq is the normal 
volume for [Z + N) nucleons, / is then varied to obtain the best fit with experimental data. 
The quantity qi^int is the internal partition function of the composite. 
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qij,int = E ( 2J k + l)e { - Ek/T) + % hCon t (2.4) 

k 

Where the summation on the right hand side is the contribution from the discrete spec- 
trum(The cut-off E max is simply the highest energy level that has been resolved for the 
given nucleus and is available from data tables); and qi j )Cor u is the contribution from the 
continuum. Without loss of generality we can write 

q AM t = J PA(E)e-? E dE (2.5) 

where we have used the abbreviation A — i + j, to stand for both i and j; Pa(E) is usually 
partly discrete and partly continuous. 

We will need both g^.mt an d Pa(E). Volumes of work are available on Pa(E). This is 
dealt with in detail in appendix 2B of |12|]. The saddle-point approximation for the density 



of states assuming a Fermi-gas model is (see eq. 2B-14 in ||12" 



p A {E) = p° A (E) x exp(/n V - a A + foE) (2.6) 



For explanations of how «o and Po are to be chosen see |T2[ . In the Fermi-gas model the 
quantity which is exponentiated is simply the total entropy S = As. Thus the density of 
states is given by a familiar expression pa(E) = p A (E) exp(S') where p\(E) is the pre-factor. 
Approximate values of p\(E) are known provided one does not have to concern with very 
low value of E (which we do need). At temperatures we will be concerned with, exp(S') in 
the Fermi-gas model is given quite accurately by exp[ir(^-) 1 / 2 }. 

In the bulk of this paper we adopt this prescription. For upto 20 F we write the density 
of state as Pa(E) = p° A x exp(S'), where the low temperature Fermi-gas expression for 5* 
as written above is used. The energy independent value of the pre-factor is fixed from 
experimentally known levels: 

Yf (2J K + l) e - E ^ T = p° A I*™* e^ E ^dE (2.7) 

fc=0 Jo 

While objections can be raised against this procedure, it achieves three objectives which 
we wanted to have: (a) we did not want to lose all information of the experimentally mea- 
sured discrete excited states ; (b) we did want to take into account the contribution from 



the continuum and (c) with this procedure calculations are fairly simple. Although, we will 
not report on all other formulae for density of states that we also used, our final results for 
the isotope populations are quite stable within reasonable variations that were tried. 
We estimate the continuum contribution as a similar integral from E max to infinity. 

qi,j,cont = / p° A p A (E)e- l3E dE (2.8) 

This process is continued upto 20 F wherein we can read off energy levels from data tables. 
For elements above 20 F, a parametrised version was used, which is given as 

,'2 /„• „-\2 \ 



Qi,j,int — exp 



(2.9) 



where W = 15.SMeV, a = 18.0MeV, k = 0.72MeV, s = 23.5MeV and e = W.OMeV. 
The first four terms in the right hand side of equation (|2.9| ) arise from a parametrised version 
of the binding energy of the ground state. The last term arises from an approximation to 
the Fermi-Gas formula for level density. This was also used in Q . For protons and neutrons 
q is 1. 

The average number of particles of a composite is given by 

{/',,) -,, Q "'- X J (2-10) 
vz,v 

However, this population is partly over particle stable states and partly over particle 
unstable states which will decay into other nuclei before reaching the detectors. 



III. SECONDARY DECAY. 

In keeping with the way experimental data are presented, we will compute ratios of yields 
of different isotopes of Boron, Carbon, and Nitrogen. To lowest order one can consider 
the (riij) obtained from equation (|2.10| ) above, remove the particle unstable fractions, and 



compare them directly with experiment. This is shown in the figures as the dotted line with 
a filled triangle plotting symbol. These populations contain only particle stable states. 



Next we consider decay of the particle unstable states. We restrict the secondary decay 
to be due to emission of six species: neutron, proton, deuteron, 3 He, triton, and alpha 
particles. Any given nucleus from a particle unstable state can in principle go to at 
most six other nuclei. As the populations are canonically distributed among the various 
energy levels, we can calculate the fraction that are in particle stable or unstable states. If 
the fraction of nuclei at the first stage in unstable states is /?-, then the number of 
nuclei left in particle stable states at the stage we call 'upto single decay' is given by 

Kj) 1 = (1 " &X"U> + £(! - flf-^fla, ]+b ^ + a, j+b ) (3.1) 

a,b l T 

where ■ is the fraction of the once decayed nuclei in unstable states. We will indicate how 
to calculate - in the next section. F a b is the width for emission by (a, b) from (i + a,j + b) 
and T T is the total width. 

We can then take these revised populations (riij) 1 and again compute the ratios. We 
label these 'upto single decay'. These are reported in the plots as the small dashed line with 
the diamond plotting symbol. Note: this is just the stable fraction of the population after 
one stage of decay, the actual population is possibly greater. 

After the first decay there may be some fraction in particle unstable states. These can 
decay, thereby, changing the population of to (riij) 2 . If we take the ratios now we get 
what we call 'upto double decay', this is denoted by the dot— dashed line and the square 
plotting symbol. Again at this stage the (n i: j) 2 represent only the sum of the stable fractions 
of the populations obtained from the initial distribution, single and double decays. 

It is clear the procedure can be continued. The fraction remaining in particle unstable 
states will continue to decrease. We found no significant difference between the 'upto triple 
decay' and the 'upto quadruple decay' calculation. Thus we do not continue beyond. Once 
again it should be noted that all the plotted populations, (n it j), (riij) 1 , (n it j) 2 , (n it j) 3 etc., 
quote only the stable fractions at freezeout, after single, double, and triple decay respectively. 

The formalism for the decay calculation is given in the next section, there the quantities 
fi,j, Ta,6 will be calculated in somewhat greater detail. The reader who is only interested in 
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the final results could jump to sections VI and VII. 



IV. THE DECAY FORMALISM. 

As the heated clusters stream out from the hot source, many of them will be in particle 
unstable states, these will decay by particle emission, for example, by emitting a neutron, 
proton, a particle etc. They will then leave a residue nucleus which may be particle stable 
or unstable; if it is unstable then it will decay further into another isotope and this process 
will continue till the residue is produced in a particle stable state. 

The primary calculation assumes that thermal equilibrium is achieved at freezeout; if 
this is true then the number of composites with % protons and j neutrons with an energy in 
the interval E and E + dE is given by the canonical factor 

dn A (E) = C hJ p A (E)e^ E dE (4.1) 

Where we have abbreviated A to mean and pa(E), from section II., is given as Pa{E) = 
p° A exp(S). The multiplicative constant p° A will, henceforth, be absorbed into the overall 
normalization constant C it j. Thus from now on the density function is given simply as 

Cij is a normalization constant such that 

PCX) 

/ C i jp A (E)e- pE dE=(n (i j ) ) (4.3) 

J 

Now of the various levels in a particular nucleus, some will be at a very low energy and 
as a result will be stable to any form of particle decay. Those that lie above an energy 

E x , v = (M x , y + Mi- x j-y - Mij) + V x , y (4.4) 

will in general be unstable to decay via emission of a particle (x, y) (i.e., a particle with 
neutron number y and proton number x) , where M X;V is the mass of the particle, M it j is the 
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p A (E) = exp 



mass of the decaying nucleus Mi- X j- V is the mass of the residue left over after decay 

and V Xt y is the coulomb barrier for that particle. Note that (x, y) could represent a variety of 
particles; in this note we will consider 'six' such particles, as mentioned in the introduction. 

As is evident from equation fl4.4|) , different particle decays have different energy thresh- 
olds. Consider an isotope as an example let us take 12 C (i = 6,j = 6). As we start 
from the ground state level and move upwards, we will encounter different thresholds. The 
lowest will be the A He decay threshold at an energy Li = E 2i2 ( in 12 C it is at 9.6MeV 
approximately) , the next higher threshold is for proton decay at L 2 = -E^o (in l2 C it is at 
18.14MeV approximately), and so on; we will get different thresholds one after the other( 
note: the order of different thresholds is different for different isotopes ). 

All nuclei of type which are formed between the ground state and the lowest thresh- 
old Li, will remain as isotopes this number is given by 

71^(0 L x ) = [ Ll C i:j p(E)e~ f3E dE (4.5) 

J 

Those that are formed between L\ and the next threshold L 2 , will all completely decay 
by 4 He emission, and these nuclei will then appear as nuclei of type (i — 2, j — 2) and must 
be added on to the population of isotope (i — 2, j — 2). 

Then, those nuclei of type which are formed between L 2 and the next threshold 
L3, will decay both by 4 He emission and by proton emission. In the next zone there will 
be three kinds of decay, and so on. We now ask, how many of the initial nuclei formed in a 
particular zone will decay by each of the channels that are available, and how many of the 
residues formed will be stable or unstable? 

To answer the above questions: we start by writing down the number of particles of type 
(x, y) with energy between (e, e + de) that are emitted, in a time interval between t and 
t + dt, by nuclei of type (i, j), lying between an energy (E, E + dE), leaving behind a residue 
nucleus (i — x,j—y) ( we may alternatively refer to (x, y), (i, j), and (i — x,j — y) by simply 
their mass numbers a, A and B where a = x + y, A = i + j and B = i — x + j — y ) 

d 3 N a = W(E } e)dedtdN(E, t) (4.6) 
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where dN(E, t) is the number of nuclei of type (z, j) initially formed at an energy (E, E + dE) 
which are still left undecayed after a time t, given by 

dN(E, t) = C i ^e- TTiE)t p A (E)e- (3E dE (4.7) 

and W(E,e)de is the Weisskopf decay probability per unit time || given by the expression 

W(E, e)de = g ala ea [a+B ^ A] Pb( e ~ B a -e) P ° B (4g) 

Pa\&)Pa 

In equation (|4.7|) , Tt(E) is the the total decay probability per unit time from an energy 
level E of the isotope A. In equation (|4.8|) , g a is the spin degeneracy factor of the emitted 
particle, 7 a is a constant of a particular decay |7]||, given by 

where m v is the mass of a nucleon. In equation ( |4.8|) , B a = Mi_ X j_ y + M ia — M it j, is the 
separation energy of the decay; <Jia+B^A} is the crossection for the reverse reaction to occur, 
( i.e. a + B — > A). It is given semiclassically for uncharged particles as 

(T[a + B^A] = nRl (4.10) 

and for charged particles as 

v [a+ B^A] = irR 2 J-^-B{e - V) (4.11) 

where R a is the radius associated with the geometrical crossection of the formation of A 
from B and a. Following the prescription of Friedmann and Lynch J7j, R a is given by 

{{A-a) 1 /* + (a) 1 / 3 ] ro , for a > 2 
ro(A-l) 1 / 3 ,fora=l 

where ro = 1.2 fm. In equation (|4.11|) , V is the Coulomb barrier for the formation of A from 
B and a. Again following fl7|, this is written in the touching sphere approximation as 



R a 



Va 



x(i-x> 2 for « > 2 

((A-o)V3 + ( )l/3) r<! > LUL « ^ Z > 

£$7S , for protons, 



where r c = 1.44/m. Also in equation (|4.8|) , Pa{E) , Pb(E — B a — e) are the respective 
density of states of the two nuclei. They have the same form as in equation ( |4.2| ). Also p° A 
and p° B are the respective multiplicative constants for the density of states, as mentioned in 
section(II). 

We note that d 3 N a in equation ( f4.6p is also equal to the number of nuclei that were 
initially formed as nuclei of type A at an energy between E and E + dE, and then decayed 
into nuclei B with an excitation energy of E — B a — e. To get the total number of states 
that decayed from a level E by emission of a particle of any allowed energy, we integrate 
over e from its minimum value V a to its maximum value E — B a , and get 

d 2 N = T a (E)C itj e- TT{E)t p A {E)e- f3E dEdt (4.12) 

where 

T a (E)= f E ' Ba W(E,e)de (4.13) 

JVa 



which on integration gives 

r.(JB)- 2lM 



Pa(E) P ° a 



( P CB p CA ^ 

(E — B a — V a ) [ —(8- 1/C) - ~^{A - l/C) 



? CB ( b3 3B 2 , 6B 6 



c v c c 2 c 3 ; 

3A 2 6A 6 



e 



CA 



3 



c \ c c 2 c 3 



(4.14) 



the derivation of the above equation is given in the appendix ), where C = n y 



i+j-x-y \l/2 



B = \/E — B a — V a , A = 0. In the above equation g a (equation (|4.8Q ), and some of the 
factors of o (equation( |4.11|) ) have been absorbed into 7^ thus 

U = la9aK&* (4-15) 

We may now integrate out the time to get 

dN a = ^^C hjPA (E)e- pE dE (4.16) 
1 T {h) 
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To get the total number of states that have decayed from nuclei of type A by channel a we 
must integrate over E from L\ to oo, 

N a = £ ^^ ) C i , j p A (E)e^ E dE (4.17) 

This integration is quite involved for as we crossover from one zone of decay (L 1; L 2 ) to 
another zone (L 2 ,L 3 ), T T (E) changes discontinously as a new channel of decay becomes 
accessible to the nuclei. Thus we break up the integration into 6 zones, corresponding to the 
6 real decay zones, and integrate within each zone independently. Note that the last zone 
extends from L 6 to L 7 = oo, and is thus considerably larger than the other zones. However, 
at the low temperatures that will be encountered, this zone will be sparsely populated. Thus 
the following approximation is valid. Within each zone, with an energy from L k to L k+1 , 
the integral can be replaced by a mean value expression, 

N a (L k , L k+1 ) = Lfc ,^ a ( Lfc+1 An itj (L k , L k+1 ) (4.18) 

L k \ L T)L k+1 

Where N a (L k , L k+ i) is the mean number of nuclei of type (or A) that were initially 
formed at an energy between L k and L k+i , and decayed by the (x,y) (or a) channel. In the 
above equation 

L k Pa)L k+1 = [ Lk+1 T a (E)C l , jPA (E)e^ E dE (4.19) 

J L k 

Of course, the left hand side is zero if channel a is not open in the region L k to L k+1 . The 
mean decay rate over all channels is 

L k (r T )L k+1 = f k+1 T T (E)C hjPA (E)e^ E dE (4.20) 
JL k 

and 

Am d (L k , L k+1 ) = [ Lk+1 C h3 p A (E)e~ /3E dE (4.21) 

Thus by summing up all the contributions from the six different zones, we get the total 
number of nuclei that have decayed from isotope A by the a channel as 
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6 

N a = ^N a (L k ,L k+1 ) (4.22) 

k=l 

To find out how many of these have decayed to stable isotopes, we must first calculate 
from equation (|4.6|) the stable decay rate T^(E). Two cases emerge in this calculation. If 
E — B a — V a > E A _ a , T s a (E) is obtained by integrating over e, from (E — B a — E s A _ a ) to 
its maximum value (E — B a ), where E A _ a is the stable level or the lowest threshold L\ of 
the residue nucleus B above which B is unstable. The expression for T s a (E) is obtained 
from that of T a (E) in equation (18) by replacing B = yj E A _ a . If E — B a — V a < E A _ a , 
then T^(E) = T a (E). Then, following a similar procedure as above for T a , we get the total 
number of nuclei A (or lying in an energy range between (L&, Lk+i), that decay by 

the a channel to a stable state as 

N*(L k , L k+1 ) = Lk fi°\ Lk+1 An id (L k , L k+1 ) (4.23) 

L k [ L T)L k+1 

The unstable decay rate from a particular level or zone is the probability of a decay per 
unit time from A to an unstable level or levels of B from which further decay can take place. 
It is easy to see that they are given simply as the difference of the total decay rate and the 
stable decay rate i.e., 

K = Ta - r* (4.24) 

The derivations and expressions for the full decay rates are given in the appendix. 

After a decay has taken place ( A — > B + a), we ask what is the population distribution 
of the residue as a function of its energy ( x = E — B a — e ). This can, in principle, be 
calculated from equation (|4.6| ) by integrating over E and e, such that (x = E — B a — e), 
the energy of the residue, is a constant. First we make a change of variables from (E,e) to 
(E, x) and then integrate over E only. We get 

dN a (x) = ( r d El ' a E ~ Ba ~ X Pfl (4g,/^) dx (4.25) 
\JB a +Va+x T T (E) p A J 

This integration is quite involved. We assume that the residue population is canonically 

distributed, but with a new temperature 1/(3' i.e., 
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dN a (x) = V id _ kj ip B (x)e 



-(3'x 



(4.26) 



There are two unknowns in this formula, the new temperature 1/(3' and the overall nor- 
malization constant T>i t j^k,i- To find these two constants we will impose that the total 
population of this interim stage (i.e. Nij^k,i )> an d the mean energy of the distribution (x), 
be reproduced by this new temperature. 

We can obtain formal expressions for the total population of the residue B as contributed 



by the decay of A, as well as its mean energy (x), from equation( [4.26D as 



N a (P',V,C) = / dN a (x) 
Jo 



i,j—>k,l 



1 + C ^* 2^ 



(4.27) 



(x(P',C)) 



1 



N a (P',V,C) Jo 



xdN a (x) 



*Di,j-*k,i 

C 2 C-K/k 



,C 2 /il3' 



+ 



+ 



1 + erf 



C 



D C 2 /Ad' 



(4.28) 



4/3 /3 1 g/3 ,7/2 y ■ 

Where the formal expression for N a (f3' ,T>,C) is used in equation (|4.28| ). The numerical 
value of N a is taken from equation( [4.22| ). The numerical value of (x) is found by explicit 
use of equation ( |4.25| ). From these two equations we obtain the two constants U it j^ k j and 
P>. 

The numerical value of (x) is derived from equation ( |4.25| ) as follows. 
1 x 



<*> = ^/ dxx(f EmaX d El > E Ba J a X PB {x) P l Ch3 e-^ 
N a Jo \JB a +Va+x lt\E) p A , 



(4.29) 



In the above equation, the numerical value of N a is taken from equation( f4.22|) . We may now 
change the order of integration to get 



i r oo r E-B a -Va E — F> — V — r n° 

The x integration is now done simply to obtain 



(4.30) 
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(x) 



N a J Ba+ v a Tt{E) 
where 1(E) is given by 

1 p% \A{E -B a - Va ) 2 e c ^ E - B «- v « 28(E -B a - y a fl^ e ^E-B a ~v, 



(4.31) 



1(E) 



Pa(E) p\ 



C 2 



C 3 



+ 



108(£ -B a - V a )e c ^ E - B *- v « 240(£ -B a - V a f/ 2 e c ^ E - B ^ 



C 4 



C 5 



^ 24Q e cVE-B a -v a 12 (E - B a - K) 240' 



(4.32) 



C 6 C 4 C 6 

In the ensuing integration over E, we, once again, replace the integral with its mean 
value expression. 



where 



^ ~WL „ dhla IT (WW C iJPA{E)e 



(4.33) 



and 



Ba+v a (rr(-E))oo = L k ( r T(E)) Lk+1 

E k >B a +V a 



Ba + V G 



'ME)), 



Ba + Va 



dEI(E)C hjPA (E)e 



-PE 



r -pG a PB_ 

Pa 



1 M 3CVnf 



/3 4 4/3 9 / 2 
C 3 0F ( / C 

F 4/3 5 + 8/? 11 /2\ 1 + er/ l v 2 v ^ 



2^ 



3 C 2 /4/3 



(4.34) 



(4.35) 



Thus the formal expressions for N a ((3',T>,C) ( equation ([4.27|) ), and (x(/3',C)) (equa- 
tion( |4.28|) ), are compared to the actual values obtained for N a (equation( fl22|) ), and (x) 
(equation( [4.33| )), and the two unknowns of equation( |4.26| ) are evaluated. We can now pro- 
ceed with further decays following the same procedure as before with decay occurring from 
a canonically distributed population at a temperature 1/(3' . 

We can thus model an n— step decay process by assuming that at each intermediate stage 
the population is canonically distributed with a new temperature and overall normalization 
constant. The decay rates to the next stage are calculated with the new temperature. 
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Following this, the fraction of the population that decays through a particular channel, and 
the mean energy of the resultant residue nucleus, are calculated. These are then used to 
secure the temperature and normalization constant of the next stage of decay. This process 
will continue till the fraction of decay to particle unstable states becomes negligible. 



V. THE CALCULATION. 

From the primary calculation we obtain that {finj) nuclei of type ( or A) are formed 
from the initial multifragmentation. The population (n»j) is distributed canonically among 
the various energy levels as demonstrated by equation( [4.1[ ). If a particular nucleus is at 
a sufficiently excited state then it will emit a particle (x, y) (or a) and leave a residue 
(i — x = k, j — y = I) (or B), which may again decay by emitting a particle (u,v) (or b) 
leaving a nucleus (k — u = m, I — v = n) ( or D ), and so on untill it finally reaches a 
nucleus (p, q) (or Z) in a stable state. We ask the question that if (riij) nuclei of type A 
were intially formed, then how many of these will finally end up as stable nuclei of type A, 
B, D ... Z. The contribution of (riij) to the final stable population of A is given simply by 
equation (|4.5| ) as 

n{ = f Ll C iJ p A (E)e- l3E dE = An^O, L x ) (5.1) 

J 

the number of nuclei initially formed as which decay to (i — x,j — y) = (k, I) is given 



as 



n A ^B = E Lk ]l a \ Lk+1 An itj (L k , Lfc+i) (5.2) 

fc=l L k \^ T/L k+ i 

The mean energy of the newly formed residue nucleus is given by 

W = la /F (TA\ An *j( B a + V a, OO) (5.3) 

We assume that this population is canonically distributed from an excitation energy of 
E = to oo with a new temperature !//?'( equation f4.26p. Extraction of the new tempera- 



ture 1/(5' and the overall normalization constant T>i^^k,i is done as detailed in section(IV). 
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In most cases, where this procedure was implemented, we obtained a new temperature 1/(3' 
which was lower than the initial temperature 1//?; however, in about 3% of the cases 1/13' 
turned out to be higher than 1/(3; this occurs when the residue of the decay process is far 
from the valley of stability. We can then calculate the number of nuclei that initially started 
out as A's and finally ended up as 'stable' -B's as 

n f A _ B = f V id ^ Kl p B {x)e- p ' x dx (5.4) 

note that in the above equation Pb{x) and Li are the density of states and lowest decay 
threshold for the nucleus of type (i — x = k, j — y = I) . This number can also be calculated 
directly by using the stable decay rates (equation( [4.23D ) 

n f M = E Lk P s Lk+1 ^(L k , L k+l ) (5.5) 

k=l Lk T/L k+ i 

The second equation is more correct as it does not depend on the assumption that the residue 
is canonically distributed. A comparison of the n A _+ B obtained from the above two equations 
gives an estimate of the error involved in the assumption of a canonically distributed residue 
population. Now we ask, what is the number of nuclei of the B's just formed which will 
decay by emitting a particle b to a nucleus of type D; this is calculated simply as 

ua^d = E Lk ^ b)lk+1 An A ^ B (L k , L k+l ) (5.6) 

fe=l L k \t T/L k+1 

the decay rates in the above equation are calculated with the temperature 1/(3'. We then 
calculate the mean energy (y) of the new distribution as 

(V) = ~ J — y» ™£ { fb An A ^ B (B b + V b , oo) (5.7) 

Using these, we continue the process on, by again calculating the temperature and norm 
of a canonical distribution, which when summed from excitation energy to oo is equal to 
n A ^B^Di and whose mean energy is equal to (y). We can then proceed to find how many of 
these will be in stable states, how many will decay on further etc. We continue this process 
till the contribution from this decay chain, A — > B — ► D — > will give numbers of nuclei 
negligible compared to the already present number in stable states. 
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VI. RESULTS OF THE CALCULATION. 



Our objective is to calculate the yields of the Boron, Carbon and Nitrogen isotopes mea- 
sured in the S + Ag Heavy- Ion collision at an energy of 22.3AMeV 0. In figures 1 to 5 the 
data are shown as empty squares. The method of calculation is simple, first we calculate 
the primary populations of the isotopes using equation ( [2 .1(J| ). We then remove the unstable 
fraction of the population, and quote only the stable part. This is denoted by the dotted 
line and triangle plotting symbol. We then incorporate secondary decay by adding on all 
the populations of nuclei that can reach a stable level of the isotopes by emitting only one 
of the six particles considered. We call these the 'upto single decay' populations and denote 
them by the small dashed line and diamond plotting symbol. We then add on all those 
unstable nuclei which can reach a stable level of the given isotopes by sequentially emitting 
any two particles of the six considered. We call these the 'upto double decay' populations 
and denote them by the dot— dashed line and square plotting symbol. We then add on all 
those that can reach the isotopes by three particle emissions, called the 'upto triple decay' 
population and denoted by the large dashed line and star plotting symbol. And finally we 
add on the 'upto quadruple decay' population denoted by the solid line and circle plotting 
symbol. As there is negligible difference between 'upto triple decay' and 'upto quadruple 
decay' we stop after 'quadruple decay'. 

To fit with experimental data, we have four parameters to tune, the obvious ones being 
the initial temperature (3 or T, the free volume Vf of the primary calculation, the ratio A/Z 
( as one does not know how much loss due to pre— equilibrium emission has taken place ) 
and an overall multiplicative constant 7i ( as we do not know how many nuclei collided in 
the experiment ). The plots are noted to be most sensitive to j3 and A/Z. Thus in fitting 
the data we first set particular values of j3 and A/Z, and calculate the multiplicities at all 
stages of decay ( Vf is varied to get the best possible fit at this temperature and A/Z ) . We 
then multiply all the multiplicities by an appropriate Ti and take the logarithm. These are 
then plotted and compared with log(counts) obtained from the experiment. We then vary 
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(3, A/Z and repeat the above procedure till a good fit is obtained. We present fits for three 
different temperatures, and different A/Z for each temperature. Vf and 7i are set to obtain 
the best fit possible for a given (3 and A/Z. 

We note that the S + Ag system is one with A = 139 and Z = 63 thus A/Z = 2.2. 
The authors of [IJj state that some pre-equilibrium emission may have taken place. As we 
do not know what proportion of neutrons and protons are lost in such a process, we start 
the calculation with the same A/Z as the S + Ag system. We start with a Z = 50 and 
A/Z = 2.2 i.e. A = 110. We start the calculation with a low temperature of 3MeV in 
Figure(|l|) ( Vf and Ti are varied to get the best fit ). We note that overall there is a slight 
excess of the heavier Nirogen isotopes as compared to data and a deficit of the lighter Boron 
isotopes, this implies that the temperature is too low and enough of the light isotopes are 
not being formed. We proceed by raising the temperature to 5MeV, maintaining the same 
A/Z. By now varying Vf and TC we find an excellent fit with the data (figure(0)). 

One may ask at this point, if there is more than one set of parameters which fits the 
data well. To answer this question we increase, first, the temperature to 7MeV, maintain 
the same A/Z and redo the calculation. We get a bad fit (figure^)). There is an over all 
deficit in the Nitrogen population and an excess in the boron population. Also we note that 
within a particular Z there is a deficit in the neutron rich isotopes. We try to remedy this 
situation by increasing the A/Z ratio. The best fit at this temperature is obtained at an 
A/Z = 2.3 (figure^)), but we still obtain an overal deficit in the Nitrogen population; the 
Carbon fit is good, but there still remains an excess in the boron isotopes especially in the 
neutron rich isotopes. 

On inspection of the fits (figures ([D to (f|)), we note that the best fit is obtained at 
figured). In this fit T = 5.0MeV, A/Z = 2.2, V fr /V = 5.5, and log(W) = 6.34. In this 
figure we note that, for the Boron populations we get an excellent agreement with the data. 
In this case there seems to be little change after single decay. For the Carbon isotopes the 
agreement is good. For Nitrogen, we have a good fit except for the case of 13 N. 

Another property of the fits noticed is that they do not seem to depend on A and Z 
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independently but rather on the ratio A/Z. As a demonstration of this, we plot in figure(j5]) 
a fit for A = 140 and Z = 63 ( i.e., A/Z = 2.22 ). We note that we are able to obtain a fit 
very similar to fig.0, with the same temperature and Vf r /Vo as in fig.|2|, but with a slightly 
lower Ti. This is very much expected, as in this case each source has a larger number of 
nucleons than before. 

VII. DISCUSSONS AND CONCLUSION . 

In this note we have presented a secondary decay formalism and performed calculations 
to fit the populations of various isotopes measured in |TJ. We obtain very good fits (fig.0) 
with experiment for the Boron and Carbon isotopes. In the Nitrogen isotopes, we obtain a 
good fit except for the case of 13 A. No particular reason could be found for this, but let us 
go over several approximations ( introduced to keep the calculation at a resonably simple 
level ) which may have contributed. 

Actual energy levels from data tables were used only upto A = 20 (equation ( |2.4|) ) for 
the primary populations. For higher masses, the emperical mass formula (equation (|2.9| )) 
was used. The secondary decay is very approximate, instead of calculating decay level to 
level, we have blurred out such details by using a smoothed level density. 

For the capture cros-section (equation( |4.11| )), we have used a simple semiclassical for- 
mula, assuming that all nuclei are spherically symmetric which is definitely not true. A 
more precise calculation involving level to level decay would use a more accurate expression 
for the cros-sections e.g., the Hauser— Feshbach formalism [H| [|TT|] . 

Still another problem lies in the assumption made in calculating the effects of higher 
order decay, that the interim populations can be taken to be canonically distributed. This 
is true only in first order decay, thus making the higher order contributions subject to some 
error. 

There is also an experimental problem according to the authors of |l[], the angular distri- 
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butions were forward peaked, indicating significant emission prior to attainment of thermo- 
dynamic equilibrium. Such an emission could affect the populations of the various isotopes. 

No doubt, incorporating changes to correct the above mentioned problems will improove 
the accuracy of the calculation. However, such changes may make the expressions analyt- 
ically intractible and one would have to resort to numerical means. This may slow down 
the calculation considerably. The calculations presented in this note take minimal computer 
time. Inspite of the shortcomings of the calculation presented above, this still remains a good 
test of the two component statistical model, and shows that such a model can definintely 
be used to explain certain experimental data quite accurately. 
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IX. APPENDIX. 



(i) Derivation of equation |4.14j . 

The full decay rate from a particular energy level E of a nucleus A (or (i,j)) which is 



decaying by emitting a particle a (or (x, y)), is given as [4.13 
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IVa 



(9.1) 



where W(E, e) is the Weisskopf decay probability per unit time given by equation (f4.8|) . On 
writing down the full expression for W we get 



E ~ B « iaP°B 



{e-V a ) exp 



7T 



[i + j-x- y){E - B a - e) 



de (9.2) 



'v a Pa(E)p° a 

Now we substitute z = (y/E — B a — e) and integrate over z and let C 
(^Vi+J -x-yj, then 



r.(s) 



Pa(E) P \ 

on carrying out this simple integration we get 



q-2z(E -B a -V a - z 2 ) exp(Cz)dz 



(9.3) 



Pa(E) Pa 



( P CB p CA ^ 

(E — B a — V a ) —{B - l/C) - —{A- 1/C) 



+ 



CB 



C 

IT 



v C 

3B 2 6B 6 N 
3A 2 6A 6 



C 



c c 2 c 3 



(9.4) 



with B = y/E — B a — V a , A = 0. The stable decay rate, i.e., the decay rate from an energy 
level E of the nucleus A to any of the allowed stable levels of B is given simply from the 
above expression by replacing the upper limit to B — JE S B where E S B is the stable threshold 
of the residue nucleus B . However in the event that E — B a — V a < E S B then the above 
mentioned replacement should not be made. In this case the total decay rate is the same as 
the stable decay rate. 

(ii) Derivation of the generic expression for L fc (r a )L fe+1 

From equation ( [4. 19| ) we obtain the definition of L fe (r a ) ife+1 as 



. (r a )i fe+1 



T a (E)C hJ p(E)e^ E dE 



now taking the expression of T a (E) from equation( |9.4| ) and substituting 



(9.5) 
B = 



{y/E -B a - V a ) we get 
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jfc \-"- a/ 



y/L k -Ba-Va 

z 3 6z 



P° A [\C 2 C 3 C 4 



-/3(z 2 +B a +V a )^ z 



now we may separate the integration into three parts 
where, 



Pa 



(9.6) 



(9.7) 



(9.8) 



L k ,L k+1 



(9.9) 



L k ,L 



'\/L k — B a — V a 

The three integrals can be done simply to give 



^ L k ,L k +\ 



C 4 /3 



(9.10) 



(9.11) 



L k ,L 



k^k+l 



C 2 (3 



L k e~P Lk - L k+1 e~^ + {1/(3 - B a - V a } (f^ - e^*" 1 )] (9.12) 



L k ,L 



k^k+l 



C 2 [3 ' 



-/3(s Q +y a ) e ^ 



+ e -PMl _ e -f3MZ +l + 



3C 2 
2/3 



/3C 2 3 3 
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} 



4/3 2 (3 C 2 J 
+ {er/(v^M fe+1 ) - er/(V^M fc )} 



(9.13) 



where M k = ^/L k -B a -V a - C/(2/3) and M k+1 = ^L k+1 - B~ a - l/ a - (7/(2/3) and C is 
the same as in equation |9T^. 
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The calculation of the stable decay rate is a bit more involved in the limits of integration 
and three cases emerge. If L k — B a — V a < E B , and L k+1 — B a — V a < E B , then 

L k (r S a )L k+1 =L k (Ta)L k+1 (9.14) 

if L k - B a - V a < E S B , but L k+1 - B a - V a > E S B , then the calculation of L k (T 8 a ) Lk+1 has to 
be done in two parts 

L k (K)L k+1 = ^ + ^)f (9.15) 

Pa 



where 



I{ = / T a (E)C hJ p(E)e-P E dE (9.16) 



the expression for this is the same as equation (|9.6|) with the appropriate change of limits. 

7| = f Lk+1 r a (E)C hJ p(E)e^ E dE (9.17) 

JE S B +B a +V a 

if however, L k - B a - V a > E S B , and L k+1 - B a - V a > E S B , then 

L k (K) Lk+1 = f Lk+1 Vl{E)C iijP {E)e^ E dE (9.18) 

J L k 

the above two integrals are rather trivial and thus detailed expressions are not presented. 
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FIGURES 

FIG. 1. Log(counts) vs Neutron number (N) — Proton Number (Z) for the three cases of 
Boron, Carbon and Nitrogen. The experimental data are from [1] S+Ag at 22.3AMeV . The fits 
show varying stages of decay for a total Zt = 50, At = 110, T = 3.0MeV, Vf r /Vo = 3.0 and 
logH = 6.82. The empty squares are the experimental data. The dotted line with the triangle 
plotting symbol is the primary calculation. The small dashed line with diamond plotting symbol 
is the 'upto single decay' calculation. The dot-dashed line with square plotting symbol is the 'upto 
double decay' calculation. The dashed line with star plotting symbol is the 'upto triple decay' 
calculation. The solid line with circle plotting symbol is the 'upto quadruple decay' calculation. 

FIG. 2. Same as fig.l but with Z T = 50, A T = 110, T = 5.0MeV, V fr /V = 5.5 and 
logH = 6.34. The best fit with the data has been obtained with these parameters. 

FIG. 3. Same as fig.l but with Z T = 50, A T = 110, T = 7.0MeV, V fr /V = 3.0 and 
logH = 6.82. 

FIG. 4. Same as fig.l but with Z T = 50, A T = 115, T = 7.0MeV, V fr /V = 3.0 and 
logH = 6.82. 

FIG. 5. Same as fig.2 but with Z T = 63, A T = 140, T = 5.0MeV, V fr /V = 5.5 and 
logH = 6.30. 
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